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Abstract 

We have implemented the newly-introduced, coherence-based technique of x-ray near-field speckle 
(XNFS) at 8-ID-I at the Advanced Photon Source. In the near field regime of high-brilliance syn- 
chrotron x-rays scattered from a sample of interest, it turns out, that, when the scattered radiation 
and the main beam both impinge upon an x-ray area detector, the measured intensity shows low- 
contrast speckles, resulting from interference between the incident and scattered beams. We built 
a micrometer-resolution XNFS detector with a high numerical aperture microscope objective and 
demonstrate its capability for studying static structures and dynamics at longer length scales than 
traditional far field x-ray scattering techniques. Specifically, we characterized the structure and 
dynamics of dilute silica and polystyrene colloidal samples. Our study reveals certain limitations of 
the XNFS technique, which we discuss. 
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I. INTRODUCTION 



Although Small Angle X-ray Scattering (SAXS) and X-ray Photon Correlation Spec- 
troscopy (XPCS) have succeeded in exploring the structure and dynamics of many inter- 
esting systems, the length scale of the observable systems is generally limited to a range 
from several nm to 100 nm, corresponding to a wavevector range of 10"'^ to 0.1 
(an angular range of 0.1° to 10°) ||l]-§||. Special difficulties are encountered when exploring 
the lower limit of the angular range, since to isolate the weak scattering from the strong 
direct beam, it is necessary to block the direct beam and extraneous scattering from slits, 
etc., which sets a boundary for the smallest detected angle. In order to probe longer length 
scales with x-rays, typically a Bonse-Hart camera is used, which can access a wavevector 
(g) range of 10~^ to 0.1 Conventional Bonse-Hart camera is one-dimensional 

coUimated which is not suitable for anisotropic samples. At a cost of reduced scattering 



intensity, there are several papers describing a 2D-collimated Bonse-Hart camera p0|-fi2 
They demand a de-convolution procedure and are inefficient in comparison to area-detector 
based method available for larger wavevectors. And the scanning procedure makes it difficult 



for time-resolved measurements. In addition, previous works [pTSj show that ultra small- 
angle could be achieved by using a very long sample-to-detector distance or a very small 
beam stop. In this case, additional interpolations with SAXS data are required because of 
the limitation of the field of view. The recently-introduced, coherence-based technique of 
X-ray Near-Field Speckle (XNFS) technique potentially offers an improved means of charac- 



teristically large length scale structure with x-rays |[T5| . In addition, XNFS is able to extend 
x-ray measurements to wavevectors (length scales) at least an order of magnitude smaller 
(larger) than may be achieved by a Bonse-Hart camera. 

The principle of XNFS is as follows. When a coherent or partially coherent radiation 
impinges on a disordered material consisting of a number of scatterers at random locations, a 
random set of phase shifts will be induced on the scattering beam. As a consequence, a grainy 
pattern will be observed in the scattered beam a certain distance away from the material. 



This pattern is called a speckle pattern [T^. In the near field region, under conditions where 
the scattered radiation and the transmitted beam simultaneously impinge upon an x-ray 
area detector, high-quality speckles can be also observed result from coherent interference 
between the incident and scattered beams. These speckles are called x-ray near-field speckle 
(XNFS) in analogy to the NFS that was initially exploited using laser sources [|17|. 

If instead of laser, one uses a high-brilliance x-ray source, it then becomes possible to 
study dense, optically turbid and/or absorbing media, in a range of length scales where 
no other x-ray or optical methods are applicable. To date, there exists a single manuscript 
describing the extension of NFS into the x-ray regime by Cerbino et al . They showed, the 
spatial power spectrum of x-ray NFS is in principle simply and directly related to the sample's 
structure factor [>S'(g)] in the range of wavevectors from 10~^ or less, to 10~^ or larger. 
Equivalently, x-ray NFS measures the density- density correlation function [5'i(r)] from length 
scales of 6 X 10^ A or more, to 1 x 10^ A or less. In addition, the evolution of the heterodyne 
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speckle pattern in time determines the sample's intermediate scattering function [S{q,t)], 
and its spatial Fourier transform, gi (r, t) . Here, we will demonstrate the implementation of 
XNFS measurements at beamline 8-ID at the Advanced Photon Source (APS) at Argonne 
National Laboratory, and will explore and discuss the utility and drawbacks of this method 
for studies of colloidal suspensions. 

II. BASIC THEORY 

In this section, we present a derivation of what we may expect to measure in XNFS 
experiments. We envision a sample with density p{x,y,z,t) and a detector located in the 
plane z = z' , so that x' and y' specify a given detector pixel. Then, we may write for the 
amplitude scattered by volume element dxdydz at [x, y, z) to the detector pixel at (x', y', z') 
at time t: 

pikR 

das = zroaoe''''-'''^o ^^^'V'^^'^'-'^^o ^(^'^''^'^'pix, y, z, t)dxdydz 

R 

where Tq is the Thomson radius, a^e^^^'^^^o '5(2^.»/>s)f^s-*:/o P{^^y,s)ds jg ^.j^g amplitude of the inci- 
dent wave at z, and 

R= \J (x' - xf + {y' - yf + {z' - zf 

6 and (3 are the real and imaginary parts of the x-ray refractive index. 

For a sufficiently uniform sample, for which the ^r-integrals of 5 and P are independent of 
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X and y, and ignoring the phase part, we may write 



dtts = iroaoe^''^ p{x,y, z,t)dxdydz—-e~^^^ 

R 



where A is the x-ray absorption length. Eq.|T] represents a quite different regime than that 
used to interpret x-ray imaging experiments. Such imaging experiments instead rely on the 
X and y dependence of 5{x, y, s)ds and /3{x, y, s)ds to create an image of the sample. 
Henceforth, we will neglect explicit mention of absorption, but otherwise take Eq. |l], as our 
expression for the scattered amplitude. 

In the near-field regime, in which scattered x-rays interfere with the incident beam within 
the coherence area of the incident beam, we are necessarily concerned with values of x' — x 
and y' — y that are on the order of 100 /im or less, and values of z' — z that are on the order 
of several millimeters or more, so that z' — z ^ x' — x and z' — z ^ y' — y. It follows that 



and, therefore. 



das — iroaop{x, y, z, t)dxdydze^^^ (3) 

z' — z 



To determine the total amplitude scattered to {x',y',z'), it is simply necessary to integrate 
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over the volume (V) of the sample, i.e. 



tts = iroao / dxdydzp{x,y, z,t)e 
'v 



ikz 



^^ik[{x'-xf+{y'-yf]/[2{z'-z)] 



Z' — Z 



(4) 



Heterodyne near-field speckle involves interference between the scattered beam and the 
incident beam of amplitude a^/ . Thus, the intensity at time t recorded at (x', y\ z') is 



I{x',y',t) = \az'f + al,as + a^'ol + ja^l^ ~ la^/]^ + a^ias + a^/a* (5) 



where at the detector a^/ = OoC*'^'^' and we have taken ^ a^/. 
Therefore, the measured intensity is 



/(x', y\ t) = |ao|^ ( 1 + iro / dxdydz 



p{x,y,z,t)- 



ik[(.'-^f+(y'~yf]/[2{z'-z)] 



-p*{x,y,z,t)'- 



'.[{^'-^f + {y'-yf]/[2{z'-.)] 



ao\ 1 — 2ro / dxdydz 



p {x,y,z,t)- 



+p"{x,y,z,t) 



[k[ix'-xf+{y'^yf]/mz'^z)]) 



cos{k[{x'-xf+{y'-yf]/[2{z'^z)]) 



(6) 



where p' and p" are the real and imaginary (absorptive) parts of the electron density, respec- 
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tively. 

Eq. ^ implicitly assumes perfect transverse coherence. To incorporate the effect of a finite 
transverse coherence length, it is necessary to introduce a mutual coherence function: 

_^^y'_y)^ e-[(-'-)V(2€^)+(.'-.)V(25^)] (7) 

where and are the transverse coherence lengths in the x- and y-directions, respectively. 
Incorporating the effect of partial coherence, Eq. |^ becomes 



I{x',y',t) ~ |ao|^ I 1 — 2ro / dxdydze 



{x'-xf (y'-yf 
2?2 + 252 



(8) 



, sin{k[{x'-xf+{y'-yf]/2{z'-z) ) 

„ cos{kUx'-xf + iy'-yf]/2iz'-z)) 

+ P {x,y,z,t)^^ -7— ^ ^ 



The first term of Eq. ^ is constant and the second term of Eq. ^ is a convolution in terms 
of real space variables. Therefore, in Fourier space, the first term becomes a 5-function at 
the origin, while the second term becomes a product. Therefore, in terms of the Fourier 
transform variables, q and p {q ^ and p 7^ 0), in the realistic case that the 2;- variations 
in the sample density occur on length scales less than which is typically hundreds 

of micrometers or more, it may further be shown that the Fourier transformed intensity is 
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given by 



~ |ao|^2ro (9) 

, / ,2 / 2,2 „2,2 \ 

-(- -z) / 1 J-x I P }y \ 

2 V(.'-.)2 + fc2e|"^(.'-z)2 + fe2^4y ^ 

^ ((^/„^)2+fc254)l/4((^,_^)2+^2^4)l/4 

{p\q,p,z,t) sin [$] +p"(g,p,z,t) cos [$])) 



where A is the thickness of the sample, p{q,p,z,t) is mixed in real and reciprocal space 
and the phase factor $ is equal to ^ f— t^t-^t^-^^t — ^^2 ^.''oL + tan"-*- 7/-* + tan~-^^^-^^ . 
Notice that the tan~^ terms in the phase factor describe an on-axis phase jump of a focused 
beam (Gouy effect) |18[ and underline the well known fact that field correlations propagate 
as the radiation field does. So these terms account for the change in the phase in the 
interference between scattered beam and transmitted beam. 

Eq. ^ which stands as a one-dimensional convolution, may be further simplified by the 
following argument: the ^-variations in p{q,p, z, t) occur on length scales set by the sample's 
structure, namely on the order of tens of micrometers or less. On the other hand, the z- 
variations in the remainder of the integrand occur on a length scale given by k/q"^, which is 
typically many hundreds of micrometers or more. Therefore, in the integrand it is permissible 
to replace each p{q,p,z,t) by its mean value, i.e. by its zero Fourier component divided by 
the sample thickness, i.e. p{q,p,z,t) ^ p{q,p,0,t)/A. This factor, which is now independent 
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of may be then taken outside of the integral, yielding 



I{q,p,t) ~ |ao|^2ro 



(10) 



-(- -z) 

dze ^ 



{p'{q, p, 0, t) sin [$] + p"{q, p, 0, t) cos [$])) 



Changing variable s = z' — z, Eq. |T0] becomes 



J(g,p,t)^|aoP^(-p'(g,p,0,t)i 



o 73n37T 7TTT7c¥ "T tan -tt^- -j- 



fz'-A 

-p"iq,p,0,t)j- 



tan" 



-1 s 



COS 



'A 



2 \^fc254+,2-^^2^ 



(s2+A;2g4)l/4^2+fc2^4)l/4 



„2_i_i.2c4 + tan + tan 



(11) 



The sample thickness A is about the diameter of the capillary equal to 0.7 mm, which is 
much smaller than the sample-to-detector distance s ranging from 53 mm to 203 mm. Thus, 
we assume that the integrand varies negligibly within the range of z of the sample. As a 
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result, Eq. |TT] can be further simplified as: 



I{q,p,t)-\ao\^^{p'{q,p,0,t) 



2 \kH^+^^kH^+^ 



(s2+fc2C4)l/4(52+fc2^4)l/4 



sm 



2 I s2+fc254 s2+fc2e 



+ tan + tan 



+ p"(g,p,0,t)e' 



„2f2 „2,2 
2 Vfc2e4+,2+fe254+p 



(s2+fc2|4)l/4g2+fc2^4)l/4 



COS 



"s2+fc2^4 



„2ii,2t4 + tan t:72 + tan 



-1 s 



[12] 



Here the term e ^ v*^^^J+^ k^^+^J (describes the effect of the partial coherence of the 
source beam on the speckle intensity. It is a product of two Gaussians with variances 
<^qx y = \J^^^^- The full width at half maximum (FWHM) of the Gaussian distribution 
is Wq^.y = 2 a/2 ln2crq^.^. Fig. |l| plots the FWHM of coherence of the source beam versus (a) 
sample-to-detector distance s and (b) coherence length ^. From Fig. [l|(a), we observe that 
Wq decreases as s increases until it reaches a constant value at a certain distance: 

So = ke, (13) 

which is the usual near-field condition, called Fresnel condition. However, for XNFS that 
requires that the way scattered radiation falls onto the sensor duplicates the actual angular 
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FIG. 1: The FWHM Wg of the Gaussian distribution for the effect partial coherence of the source 
beam, plotted versus (a) sample-to-detector distance for different coherence length and (b) coher- 
ence length for different sample to detector distances. In (a), dashed lines denote the boundaries 
separating far- field and near-field for 8 fj,m and 163 fim from left to right. In (b), the intersections of 
dashed lines and color lines denote the differences of the aq at different s for ^ = 8 /xm and 163 /im. 
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distribution of the scattered intensity, a much stronger condition should be satisfied [|T^] 



(14) 



where a is the size of the scattering particle. With this condition, the source beam fills the 
whole field of view and the speckle size is related to the actual size of the probing material. 
In Fig. [l|(b), Wq is plotted versus coherence length for different values of s. For the estimated 
coherence lengths at SID indicated via the dashed lines in (b), in principle, one expects to 
observe an s-dependent change in the 2D speckle intensity. However, in reality, due to the 
limited spatial resolution, which in turn limits the q range, and more critical, the sensor 
response [0, we have not been able to observe this s-dependent variation, as we discuss in 
more detail below. When <^ k'^^t ^^"^ ■^^ ^ k'^^y expect at 8-ID, then, Eq. |12] may 

be simplified even further: 



2ro -—(-^+^ 
I(q,p,t) ^ \ao\'^——e ^ V'^^^^ '^^^y 
k 



(15) 



(p'(9,P,0,t) sin 
+ p"{q,p,0,t) cos 



{q^+p'^> _^ £ ( J_ _(_ J_ 



1 (_{f+f)s 

2 I k 



fc2, 



where the term e ^ K'^^^x k^iyj jg introduced as the spatial coherence transfer function in 
Ref. [0. 
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Introducing = tan ^{p'/p"), which is always small, except at x-ray energies near an 
absorption edge, and p = (p')^ + (p")^, we may re-write Eq. |T3| as 



2ro - 
/(g, p, t) ^ Io^p{q, p, 0, t)e 



2 \^ k'^^x ^'^^y 



(16) 



sm 



k k k^^ k^y 



2r 



where p and q are the wavevectors obtained in the x— and y— directions, respectively, by 
numerically Fourier transforming the CCD image, Iq = |aoP, A is the sample thickness, p 
is the electron density in Fourier space, and where the latter equality defines the transfer 
function T{q,p). It is worth emphasizing that the transfer function T{Q) with Q = {q,p) is 
written as 



T{Q) 



9 2 2 

2 V 1.2 12 ^ ,.2f2 ) ■ 



(17) 



Eq. |T6| immediately allows us to calculate the static structure factor S{Q) in terms of 
measured quantities. Specifically, 



^(40) = (|p(g,p,0,t)|'> 



I{Q,t) 



(4ro%^|T(Q)PA^) 



14 



Similarly, it is straightforward to show that, in the context of XNFS, the normalized IFS is 



9\{Q.^) = — ^T7^\ (19) 



(p(g,p, 0,t)p*(g,p,0,t + r)) 

|p(g,p,o,t)|'> 



IiQ,t)' 



According to Eq. 19, (71 is independent of the transfer function T{Q), and in turn does 
not depend on sample-to-detector distance, s. 

In summary, the intensity measured in the near field speckle experiments is proportional 
to the density of the sample rather than the modulus squared of the density as in conventional 
far field speckle experiments like XPCS and SAXS. Thus, the time autocorrelation of I{Q) 
gives directly the intermediate scattering function gi (Eq. When the delay time r is 

chosen to be zero, we obtained a quantity proportional to the static structure factor S{Q) 
(Eq. |I8]) times the NFS transfer function \T{Q)\^. 



III. DETECTOR DESIGN 

High spatial resolution and high detection efficiency are key goals for imaging the x-ray 
speckles in XNFS experiments. In XNFS experiments, in order to resolve micron-sized par- 
ticles, it is necessary to employ detectors capable of resolving micrometers. A typical x-ray 
imaging detector consists of a crystal x-ray scintillator, a microscope objective and a fast. 
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high-resolution, large-dynamic-range charge-coupled device (CCD)-based camera. The scin- 
tillator converts x-rays into visible light; the objective collects the visible light and magnifies 
the visible-light image; and finally the CCD camera records the image. Our aim was to de- 
sign a detector with the best combination of scintillator and objective to achieve the optimal 
combination of spatial resolution and detection efficiency for XNFS experiments. 

A key characteristic of an objective is its numerical aperture (NA). The numerical aperture 
(NA) of an objective defines the largest angle of light acceptance as well as the light collecting 
power. The detection efficiency scales as (NA)^. Thus, a large NA objective is necessary for 
high detection efficiency. It is common to use immersion oil of high refractive index (n = 
1.515) between the front lens of the objective and the scintillator to achieve high numerical 
aperture. We employed an Nikon Plan Fluor 40x oil immersion microscope objective with a 
numerical aperture of NA = 1.3, a working distance of WD = 0.2 mm and a field of view of 
diameter 0.67 mm. This objective uses an infinity focused optical system with a reference 
focal length of 200 mm. In our case, with the implementation of a tube lens with adjustable 
focal length from 25 mm to 150 mm, the 40x objective gives a real magnification of from 5 
to 30 times. 

Generally, taking into account both the effects of diffraction and depth of focus, the spatial 
resolution [R] function of NA is given by pO[ : 



R = [{p/NAf + (gA^NA)2]V' (20) 
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FIG. 2: Spatial resolution (/xm) versus numerical aperture NA of the objective for different x-ray 



absorption lengths of the scintillator adapted from Ref. |20|. 



where, p = 0.18 and q = 0.075 are constants obtained by numerical simulations ||20|| and A 



is the x-ray absorption length of the scintillator. Based on Eq. ^ one can plot R versus 
NA for different Az, as shown in Fig. ^ It is clear from Fig. ^ that to achieve a spatial 
resolution of a micrometer or less with high detection efficiency (NA > 1.0) , we have to 
choose a scintillator with x-ray absorption length of 10 fim or less. 

Besides the x-ray absorption length, several additional characteristics of the crystal scin- 
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TABLE I: Characteristics of the scintillators 

scintillator x- light refraolrasitjWavelength 

ray yield in- (g ■ cmol^) 

ab- per dex max. 
sorp- keV emis- 
tion sion 
length (nm) 

(aH 

LYSO 9.5 32 1.81 7.1 420 
(Lui.gYo.sSiOsiCe) 

YAG:Ce 25 8 1.82 4.55 550 
(YaAlsOisiCe) 

CdW04 6.5 12 to 2.2 7.9 475 
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tillator are critical for x-ray imaging, including: high x-ray stopping power; high light yield; 
the emission wavelength being compatible with CCD readout (400 nm - 700 nm); a similar 
refractive index to the immersion oil (n=1.515), and a small thickness to minimize spherical 
aberration. At a minimum, the scintillator thickness should be smaller than the working 
distance of the objective (0.2 mm), so that the objective can focus to the upstream side of 
the scintillator. 

As shown in Table |, one potential candidate is YAG:Ce (Y3Al50i2:Ce) which has been 
widely used in x-ray imaging detectors. However, its x-ray absorption length at 7.44 keV is 
about 25 yUm, which is not suitable for submicrometer resolution detector with a NA = 1.3 
objective (Fig. |^). Another candidate is CdW04 whose x-ray absorption length at 7.44 keV is 
only 6.5 fim. However, It is very difficult to obtain thin crystals of CdW04. It easily breaks 
before being thinned to the desired thickness (< 0.2 mm) due to its (010) cleavage plane. 
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In addition, the CdW04 refractive index (n=2.2) is much different from that of immersion 
oil (n=1.515), which will induce a relatively large spherical abberation. In our setup, we 
use LYSO (Lui.8Yo,2Si02) which appears to be the most appropriate candidate overall. Its 
x-ray absorption length at 7.44 keV is 9.5 fim, slightly larger than the x-ray absorption 
length of CdW04 but still good enough to produce high spatial resolution. Its light yield 
is 32 photons per keV, much higher than both CdW04 and YAG:Ce. Its refractive index 
(n=1.81) is not too far from that of immersion oil inducing less spherical abberation than 
CdW04. In addition, we found a manufacturer providing two-sided polished LYSOs with a 
thickness of 0.15 mm, just appropriate for our objective, although, as we show below, even 
thinner would lead to reduced spherical aberrations. 



Here, following Koch |21|, we briefly discuss spherical aberrations of our system. Gen- 
erally, "spherical aberration" defines the aberration that is introduced when the light from 
near the center of the lens has a different focal length from light from the margins of the 
lens. We begin with calculating the spherical aberration induced by a general material with 
refractive index of rii and a thickness of t. According to snell's law, an incident ray with an 
angle of 6i against the optical axis will refract at the interface of the material and another 
media with refractive index n2 with a new angle 62: 

rii sin 9i = n2 sin 62 (21) 
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Then, the distance from the virtual source of ray and the interface is given by 



rf = t^ = t!^^ (22) 
tan 6*2 ui cos 9i 



So the shift of the focus is given by the difference between d for the paraxial ray and marginal 
ray: 



Az = d{e,, 02 = 0)- d{e,, 02) = t- f 1 - ^ ) (23) 

rii \ cost/i 




where we used tan^i/ tan ^2 = d/t, t is the thickness of the scintillator, sin ^2 ~ tan 6*2 = 
a is the aperture of objective and p' is the normalized radius of the ray on the entrance pupil. 
Following the addition theorem for the primary aberrations |^ (p246): Each primary aber- 
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ration coefficient of a centered system is the sum of the corresponding aberration coefficients 
associated with the individual surfaces of the system, we can replace p' by the radius of the 
ray on the exit pupil p, and obtain the wave aberration due to shift of focus: (p550): 



*fe)4(;)^A./4(;)%'(i<^=^) (24) 



And the spherical aberration is represented by |^ (p530) 



^|J{q)s, = A'o,oP^ (25) 



Hence, from Eq E4l we find that Aq^q is of the form: 



= l^^'''^^ (26) 



where we use NA = a/r. Born and Wolf ||2^ also derive the tolerance condition of the 
spherical aberration: 

v4'o4o < 0.94A (27) 



Hence, the maximum tolerated thickness of the material satisfies 



t<0.94A/(-NA^^^^^i^^) (28) 



2 2nf 
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For rii = 1.515, n2 = 1 and A = 420 nm, the tolerated thickness is 1.5 /im, which is too 
hard to reach. So the correction of spherical aberration is required. The objective we used 
is an oil-immersion objective with a working distance of 0.2 mm that is corrected for glass 
coverslips of a thickness of 0.17 mm. The spherical aberration induced by using a different 
refractive index n' can be compensated by choosing a different thickness t', which is given 

byfH 

t = ^3 29 

When working with a crystal scintillator of a refractive index of 1.81 and a thickness of 0.15 
mm, the equivalent thickness of the glass and oil is 0.155 mm. So what we need to do is 
remove the coverslip and increase the oil thickness to 0.215 mm, then the resultant spherical 
aberration is within the tolerance. 

The lateral displacement of the image for the spherical aberration is given by ||22| (p531): 



Then the resolution due to the spherical aberration is proportional to the lateral displacement 
of the image 

R CX Ayma. = ^AYrr^a. = ^^040 ~ ^^A^ (31) 

where AYmax is obtained when Y = a, and we used sin a = a/R and the sine condition 
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y ■ nsin6' = F sina. So the contribution of the spherical aberration is not trivial for high 
numerical aperture objective. To lessen this effect, one can either choose small NA objective, 
or choose a scintillator with thinner thickness. 

IV. EXPERIMENTAL SETUP 

X-ray near-field speckle (XNFS) measurements were carried out at beamline 8-ID-I of 
the Advanced Photon Source (APS) at Argonne National Laboratory, using x-rays of energy 
7.44 keV and half the source beam size available at 8-ID-I hutch - about 0.5 x 0.5 mm^. 
Fig. ^ shows a sketch of the optical setup located in the 8-ID-I hutch. From right to left, 
we have the beam source, the sample stage, the scintillator, the microscope objective and 
the CCD camera. In XNFS experiments, no spatial and spectral filtering of the direct beam 
are required [|T^. We inherit the exiting setup for XPCS which gives an energy resolution 
of / E ^ 3 X 10~^ and remove all the slits letting the full beam impinge onto the sample 
and record the interference pattern of the transmitted and scattered beams by means of our 
detector placed at distances from the sample ranging from z = 53 mm to 203 mm. 

As described before, we use two-sided polished LYSO (Lui,8Yo.2Si02) with a thickness 
of 0.15 mm to convert x-rays into visible light and a Nikon Plan Fluor 40x oil immersion 
microscope objective with a numerical aperture of NA = 1.3 to magnify the image. Both 
of them are mounted on a piezo electric stage which has a mechanical manual adjustable 
coarse travel range of 4 mm and a piezo electric travel range of 20 /im with a resolution of 20 
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FIG. 3: A sketch of the optical setup for XNFS experiment. The key components from right to 
left are: the beam source, the sample stage, the scintillator, the microscope objective and the CCD 
camera. 

nm. As a result, by carefully tuning the distance between the objective and the scintillator, 
we are able to focus the image which is in turn magnified 5 to 30 times by a tube lens 
with adjustable focal length and then recorded by a "CoolSNAP" CCD camera. The camera 
features a 1392 x 1040 pixels with size 6.45 fim x 6.45 fim and a maximum frame rate 
of 56 Hz. All the images so-obtained are subsequently cropped to 1024 x 1024 pixels for 
convenience of 2-D Fourier transformations in the data analysis. 

There are several contributions to the detection resolution: (1) The resolution of the 



scintillator - for x-rays incident onto the scintillator at a single point. According to Refs. pO 



21| , this is typically 0.1 /im for 7-8 keV x-rays, smaller than optical limits on the resolution. 
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(2) The resolution determined by the diffraction hmit and the defect of focus of the objective 
and the scintillator via Eq. With an objective with NA = 1.3 and LYSO with an x-ray 
absorption length of 9.5 fim, we obtain a spatial resolution of 0.98 fim. (3) The reduction in 
resolution caused by the spherical aberration due to the use of s scintillator of mismatched 
refractive index, which is proportional to the thickness of the scintillator and the cube of NA 
(Eq. ^T]). Since our objective is corrected for spherical aberration and the additional spherical 
aberration induced by replacing the glass coverslip by the scintillator can be compensated 
for by adjusting the oil thickness , this factor is not critical to the spatial resolution. (4) 
The (demagnified) size of the CCD pixel, if too large, could limit the resolution. For a 
magnification of 30, the CoolSNAP can resolve a length scale as small as 6.45/30 = 0.215 fim 
which does not limit the resolution. As a result, the spatial resolution of detector should be 
largely determined by factor (2), which is about 1 fim. 

Hence, we are able to estimate the maximum q range. With r^in = 0.98 fim, in reciprocal 
space, Qmax = ^ ~ ^ X 10"^ A~^. On the other hand, the lower limit of the wavevector g^m 
should be determined by the largest accessible length scale. In principle, the largest length 
scale is the size of the measured scattering image equal to 1024 x 0.215 fim = 0.22 mm, so 
Qmin = Q 22mm ~ 3 X 10~^ A~^. Howcvcr, practically our data shows identical low q profile 
that is independent of which sample is being studied and of the sample-to-detector distance. 
This is likely due to the beam structure on long length scales. So, the realistic useful g^m is 
of the order of 10^^ A~^. Nevertheless, the g-range achieved is at least a decade below the 
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range accessible to the conventional XPCS experiments. 



V. SILICA 0.45 nm SUSPENSION 

The first sample, we measured at 8-ID-I, as an initial test of our XNFS setup is a colloidal 
suspension of silica particles of 0.45 /im diameter and a volume fraction (0) around 0.05. 
The static structure factor peak of this sample is expected to be located around 10""^ 
which means that the sample has a uniform scattering profile in the q range accessible in our 
XNFS setup {i.e. S{q) ~ const). In other words, the intensity profile /(g) we measured from 
this silica suspension should simply result from the transfer function T{q) (Eq. |16]). Thus, 
it should be an ideal sample to examine the sample-to-detector distance (s) dependence of 

nq). 

Illustrated in Fig. |^(a) is the typical raw image of the silica suspension measured at 
z = 53 mm. The image contains 1024 x 1024 pixels with a pixel size of dpi^ = 6.45 fim. The 
magnification of the detector is set to 30. Thus this image corresponds to a region of size 
0.22 mm x 0.22 mm in the sample. The speckle pattern appears quite obscure and weak 
due to the large-scale, static background. After averaging over 1000 frames, the speckle 
pattern is washed away leaving the static beam profile fluctuation unchanged, as shown in 
Fig. H(b). In order to remove the large-scale static fluctuations, we perform a normalization 
of each raw image by dividing each one with an average image, averaged over 1000 frames, 
as shown in Fig ^(b). Fig. ^(c) presents one example of the resultant image, which reveals a 
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FIG. 4: (a) A raw single frame of scattering image of Si02 with diameter of 0.45 ^m at sample- 
to-detector distance of 53 mm. (b) The image averaged over 1000 frames, (c) Normalized image 
obtained by dividing raw image with averaged image. Each image consists of 1024 x 1024 pixels, 
(d) Two-dimensional fourier transform of (c). 
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clear, uniform speckle pattern Next, a two-dimensional Discrete Fourier transform 

is performed via 

^ N-lN-l 

Hq^p) = j;^'^'^Hx,y)^wi-ii<ix + py)/N) (32) 

which produces the g-space image, as shown in Fig. §(d). Here, N=1024. As a result, the 
corresponding q coordinates are given by: 



2n X 2n y 

" P = TT^l^ (33) 



NArAr NArAr 

where Ar is the size of one pixel, equal to (ipjx/30 = 0.215 /im. 

We report in Fig. ^ several examples of the magnified Fourier transformed image corre- 
sponding to the region inside the dashed lines in Fig. ^(d). Different panels are obtained 
at different sample-to-detector distances: (a) 53 mm, (b) 103 mm, (c) 153 mm and (d) 203 
mm. In each image, prominant fringes can be seen. It is clear that the fringes become 
finer when the detector is moved away from the sample. This agrees with the theoretical 
prediction that the transfer function |T(g)p is proportional to a sine term whose frequency 
depends on the sample to detector distance (Eq. |17D. Note that the rings of the Fourier 
transformation are not azimuthally uniform. This effect may be related to asymmetry in 
the coherence of the source beam. However, the envelop of the asymmetry has no obvious 



s-dependence, in contrast to what may be expected on the basis of Eq. 12. Thus, we do 



not understand this asymmetry in detail. In fact, examination of these data (Fig. in the 
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light of Eq. |T2|, suggests that the predicted effect of a finite coherence length is not playing 
a role in determining these data, presumably because the width of the Gaussian in Eq. is 
greater than our accessible g-range of Qmax ~ 2 x 10~'^ A^^ even for the largest values of s 
studied. In Ref. the authors showed that by far the largest contribution to the q decay 
of the speckle power spectrum is due to the sensor transfer function. It is highly possible 
that this asymmetry rings are due to the sensor response. 

To quantify |T(g)p, we plot in Fig. |^ the azimuthally-averaged intensity profile (sym- 
bols) versus wavevector [q) for different values of s varying from 203 mm to 53 mm with a 
decrement of 10 mm. The thick solid lines in Fig. ^ are the fits of I{q)s to the relation: 

liQ) = A^ \TiQ)\' + ^2e-(^-«)'' + B (34) 



with Q = \Q\ = yq'^ + p"^ and T{Q) as a simplification of Eq. ITT 



T{Q) = e"('"-«) ■ e~'' i^'^. + kH^) . gin + ^ + y,)] (35) 



which consists of a product of three terms. The first factor is additional to Eq. in order 
to take into account the limited spatial resolution of our optical system. As in Eq. [T^, the 
second term derives from the partial coherence of the incident beam. How the coherence 
length enters is somewhat counterintuitive: It is proportional to the width of the Gaussian 
term in reciprocal space. However, this term is set to unity for fitting, since it plays no 
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significant role in the accessible q range, as we observed in the context of Fig. |^. The 
third term is a sine function that describes the fringes produced by the interference of the 
scattered beam and incident beam. In x-ray domain, this term is called phase contrast 



transfer function ||T5|. Note we assume ix = iy = i for the simplification of a phase factor 
in the sine function. 

Besides the contribution of T{Q) and background noise, the second term in Eq. |3^, in 
the form of a stretched exponential decay, describes a combination of contributions from the 



sensor transfer function |19| and multiple scattering. The existence of multiple scattering 
is evident based on three observations. Firstly, the sine square term goes to zero period- 
ically, which should make the intensity profile exhibit minima with the same magnitude. 
However, the measured minima decrease with q. The second piece of evidence pointing to 
the importance of multiple scattering comes from the dynamic data (see later), which dis- 
plays a g-dependent decay rate and exponent that mirrors and anti- mirrors the form of T{Q) 
(Fig. |lO| and Fig. |lT]). This indicates that we are measuring faster dynamics due to multiple 



scattering where single scattering vanishes. Thirdly, very strongly scattering samples- i.e. 
deliberately multiply scattering samples - show no minima at all. 

Hence, we fit the intensity data in two steps. In the first step, we focus on multiple scatter- 
ing. Note that the single-scattering term vanishes at gs satisfying | + j§s + 
where n = 1, 2, 3.... As a consequence, we extract the intensity at the Q values corresponding 
to those minima in I{Q) and fit them with only the multiple scattering term plus the back- 



nvr. 
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ground. The characteristic length v is fixed to 500 yum. The amphtude A2 and the exponent 
P were allowed to vary. The resultant multiple scattering term is plotted as the thin solid 
lines in Fig. ^. The best fit exponents /3 are plotted in Fig. ^ versus the sample-to-detector 
distance s. /3 fiuctuates around 0.345, giving an empirical stretched-exponential form for 
the intensity of the multiple scattering. Next, the remaining intensities after subtracting the 
filled multiple scattering and background are fitted with the transfer function (Eq. ^5]). We 
use the measured sample-to-detector distance s and try to find one set of w, and if that 
works for all the s. The only varying parameter is the overall amplitude. The set w = 1.8 /xm, 
^ = 163 fim and ip = 0.026 yields a good fit for all values of s studied, as shown by the 
thick solid curves in Fig. ^. The value of w provides an estimation of the resolution of our 
detector, and is close to but somewhat larger than our estimate 1 fim. The small value of ip 
indicates that the x-ray absorption for this silica sample is essentially small. In general, the 
fitting reproduces the data with few fitting parameters, confirming the theoretical relation 
between the sample-to-detector distance s and the transfer function T(g), and consequently 
confirming the feasibility of our experimental setup. 

With the same principle of XPCS, the fiuctuations of near-field speckles should refiect the 
dynamics of the sample. As shown in Eq. |l^, the time autocorrelation of the intensity gives 
rise to gi instead of (72 in XNFS experiments. Hence, we presented in Fig. p, the normalized 
intermediate scattering function (5'i(t)) versus delay time (r) for r between 0.4 s and 319 s 
for different s at g = 8.180 x 10"^ A~^. The gis collapse into one curve for different s, which 
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agrees with the theoretical prediction that gi has no s-dependence due to the cancelation 
of T{q). However, for a larger wavevector of g = 1.420 x 10^^ the gis do not overlap 
for different s, as shown in Fig. ^. To elucidate the reason for this discrepancy and quantify 
our observations, we have fitted gi measured at different s and g to a stretched exponential 
form: 

91 = e-(^*)" (36) 

The best-fit relaxation rate (F) (Eq. |36|) versus wavevector q is illustrated in Fig. |10]. The 
values of F at successive s are displaced by a factor of 1.1 from the previous s value for 
clarity. Generally, F(g) at different s show a behavior, illustrated by the dashed line in 



Fig. [ly. However, peaks are observed at the q positions coinciding with the q positions of 
the dips in the transfer function T{q) (Fig. |^). 

Away from these multiple scattering peaks, F(g) increases as g^ versus g, which is rea- 
sonable for a Si02 suspension undergoing Brownian motion. Quantitatively, for Brow- 
nian motion, we expect F(g) = Dmq^- As a result, we derive the value of diffusion 
coefficient D^n = T/?^ ~ 1.167 x 10~^^m^/s. According to the first order hydrody- 
namic interactions Dm = Dq{1 + 1.450) ||23|, where Stokes-Einstein diffusion coefficient 
Dq = 1.048 X lO^^^m^s^^ with the Boltzmann constant, T the room temperate 

equal to 293 K, t] the dynamic viscosity equal to 1 x 10\g/ (m ■ s), we obtain (p ^ 0.07, which 
is reasonable. 

At the g positions of the peaks, where the single scattering amplitude goes to minimum 
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because of the zeros of sine^ term in T(g), we hypothesize that we are measuring muhiple 
scattering of the sample. This theory explains why we obtain faster dynamics at those q 
positions|2^. Fig. |1T] shows the corresponding best fit exponent a, which exhibits similar 
fluctuation patterns as T(g) and supports our hypothesis. Underlying this hypothesis is the 
idea that the rapid variations of T(g) versus q may be associated with single scattering, 
whereas the intensity of multiple scattering likely shows a relatively smooth g-dependence. 
Accordingly, if, for a particular set of data, the scattering minima due to T{q) are indistinct 
(do not send the scattering intensity to zero) then it follows that the XNFS data set in 
question suffers from multiple scattering. Hence, to calculate (71, we have to pick q smaller 
than the first dip of T{q) so that the measured coherent scattering is reliable. 

To further test this idea, we carried out measurements on a sample that could be expected 
to show very strong scattering and therefore strong multiple scattering, namely a 3 mm-thick 
sample of Gillette Foamy shaving foam, which is know to consist of a dense foam of micron- 



sized air bubbles in aqueous liquid. Fig. [12] shows the scattering intensity from such a sample 
as a function of g, obtained using the XNFS prescription. However, in contrast to the more- 
weakly scattering silica spheres, discussed above, evidently, in this case there are not the 
oscillations in intensity that are expected for XNFS, i.e. there is no evidence that the XNFS 
|T(g)p is displayed in these data. We infer that this is indeed the result of multiple scattering 
and that the x-ray scattering from 3 mm-thick foamy is completely in the multiple scattering 
regime. This implies that |T(g)p is a signature of single scattering. We can also calculate 
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Qi for foamy according to the XNFS prescription. This is shown in Fig. |T3[ The dynamics 
is pretty slow, on the order of 0.1 s~^. 

These results point to another difficulty with the XNFS method (which is common to 
ultra-small-angle x-ray scattering methods in general) namely that multiple scattering must 
be carefully considered and if possible eliminated. In the case of foamy, a sufficiently thin 
sample (much thinner than 3 mm) would have eventually had reached the single scattering 
regime. Interestingly, in the case of XNFS, in contrast to more traditional USAXS meth- 
ods, the existence or not of multiple scattering may be straightforwardly and immediately 
recognized from the intensity profile i.e. |T(g)p as we discussed previously. 

VI. POLYSTYRENE 4 SUSPENSION 

In this section, we present the XNFS data obtained from a colloidal suspension of 
polystyrene particles of a diameter of 4 fim. This sample is not as stable as the last sample, 
since particles with 4 /im undergo sedimentation. The static structure factor peak of the 
polystyrene suspension of this size lies within the g-range accessible by our XNFS setup. 
Hence, we expect to observe more complicated intensity profiles with the contributions from 



both structure of the suspension and the transfer function. Illustrated in Fig. |TJ are the scat- 
tering intensities (symbols) plotted versus q obtained by azimuthally averaging the fourier 
transformed scattering images averaged over 1000 frames for sample-to-detector distances 
s = 113 mm, 143 mm, 173 mm and 203 mm (from top to bottom). The corresponding 
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transfer functions (T(g)) obtained by fitting data of silica sample measured at the same s 
are plotted in solid lines with the same colors for easy comparisons for the peak positions 
of Ttaib- The intensity data deviate from the lines. Firstly, the peak positions of the data 
matches the ones of T{q). There might be one extra peak located at q around 10~^ for 
all s, which is coming from the static structure factor peak. In addition, the peaks of this 
sample are less sharp than the ones of silica sample. We attribute this to the fact that the 
scattering contrast of polystyrene is less than that of silica, leading to weaker scattering for 
the polystyrene sample, even though the PS particle are larger. 



Illustrated in Fig. |T5] are the normalized intermediate scattering functions (gi) at q = 
8.18 X 10~^ for delay time from 0.02Doq^ to 2QDQq^ seconds and sample-to-detector 
distances s from 203 mm to 113 mm with an interval of 10 mm. The for the polystyrene 
suspension do not totally overlap for different s at this q position, but decays slightly faster 
when the detector moves closer to the sample stage. We reason that this is the result of 
the sedimentation of the polystyrene particles, which leads to denser sample with faster 
dynamics. 

Following the same procedure as for the silica sample, the best-fit relaxation rates (P), 
obtained by fitting the one of the 100-frames gi with single exponential form (Eq. 



are plotted versus wavevector q for different s in Fig. [T6|. The values of P at different 
s are displaced for clarity. Similarly, peaks that correspond to dips of transfer function 
are observed, confirming our conclusion about the measurement of multiple scattering at 
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the minimums of transfer function. In this case, the peaks are more visible than the ones 
observed in sihca sample, indicating stronger multiple scattering in this sample with bigger 
polystyrene particles. 

VII. FUTURE WORK AND CONCLUSION 

In conclusion, we have presented the implementation of the new coherent x-ray technique - 
X-ray Near Field Speckle as well as its applications and limitations. Clearly, XNFS is capable 
of obtaining ultra small angle x-ray scattering and x-ray photon correlation spectroscopy with 
its simple setup and direct relationship to the density correlation function. It effectively 
extends to wavevectors an order of magnitude smaller than the wavevector range covered by 
conventional SAXS and XPCS, and enables us to explore the static and dynamic structures 
of micrometer-sized samples. We believe this technique will be valuable for optically dense 
and turbid samples which induce strong multiple scattering optically. 

Technically, XNFS is not difficult to realize. Speckle pattern is produced simply by 
letting both scattered beam and transmitted beam impinge onto the detector. It does not 
require the spatial filtering as did in XPCS, which allows us using the whole source beam 
and in turn greatly enhance the speckle contrast. As a consequence, it does not require 
laborious alignments. All the efforts were devoted to the design of the detector. High 
numerical aperture objective was employed to produce high spacial resolution and efficient 
light collection. The measurements gives convincing results, which proves the feasibility 

36 



of this setup. Improvements could be made on several aspects. One is to utilize thinner 
scintillator, which will give rise to less spherical aberrations. A faster CCD camera will for 
sure improve the probing range of the dynamics of this technique. 

One key difficulty of this technique is due to the transfer function T{q). It entangles with 
S(q). it is straightforward to characterize the structure factor peaks and dips located smaller 
than the q position of the first dip of T{q). However, this would make the reliable q range very 
small. If XNFS is to realize its full potential, it will be necessary to figure out an effective 
way to deconvolve the static structure factor from the transfer function in the future. One 
possible way is to use as small as possible sample-to-detector distance though with a cost in 
scattering contrast. Strong absorption samples might not be affected by this factor due to the 
phase factor induced in the sine term of T(g)(Eq. ^). Another possible improvement might 
be made by measuring a control sample with exactly the same material but uniform S{q) in 
the accessible q window, then dividing the intensity profile of the interested sample and the 
control sample. Another important issue is multiple scattering. Evidence of the existence of 
multiple scattering comes from the intensity profile and sample-to-detector dependent decay 
rate. Making the sample as thin as possible should solve this problem. 

Finally, the limited dynamic range of the CCD (12 bits) means that weakly scattering 
samples cannot be studied, because they give a scattered intensity that is less than 1/4000 
of the direct beam intensity and so their scattering contribution cannot be reached. Using a 
16, 18 or 20 bit CCD would extent the range of possible XNFS sample correspondingly. 
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FIG. 6: Intensities of Si02 suspension at different sample-to-detector distances. Symbols are the 
data. Thick solid lines are the fittings based on Eq. ^ Thin solid lines are obtained by fitting the 
local minima with a stretched exponent, describing the contribution of multiple scattering. Data 
are displaced by a factor of 10 for clarity. 
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FIG. 7: The best fit exponent /3 versus sample-to-detector distances varying from 53 mm to 203 
mm for silica 0.45 fim suspension. 
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FIG. 8: Normalized intermediate scattering functions of Si02 suspension versus delay time r mea- 
sured at different sample-to-detector distances and q = 8.180 x 10~^ A. 



45 









1 - 


c 


- 


o 




o 


! 


c 


0.8 - 




Ll_ 






- 


c 




(D 






0.6 - 


D 




O 




m 


- 


(D 




D 


0.4 - 


13 


£ 








c 


0.2 - 


13 




CD 








D 




E 


- 


o 


z 





q = 1 .420e-004A 




10 



10 



10 



w 


niQ+ n^? mm 

L^lo 1 UOO 1 1 M 1 1 




niQ+ Di^'^ mm 


w 


niq+ ny^ mm 

L-.'lo 1 / \J 1 1 M 1 1 


• 


Dist 083 mm 


o 


Dist 093 mm 


o 


Dist 103 mm 


o 


Dist 113 mm 


o 


Dist 123 mm 


o 


Dist 133 mm 


o 


Dist 143 mm 


o 


Dist 153 mm 


o 


Dist 163 mm 


• 


Dist 173 mm 


• 


Dist 183 mm 


• 


Dist 193 mm 


• 


Dist 203 mm 



10 



10 



Delay times (s) 



FIG. 9: Normalized intermediate scattering functions of Si02 suspension versus delay time r mea- 
sured at different sample-to-detector distances and q = 1.420 x 10~^ A~^. 
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FIG. 10: The best-fit decay rate T versus wavevector q for silica of diameter 0.45 /xm measured at 
different s. Points are tlie fitting results and the lines are guides to the eye. Data are displaced by 
a factor of 1.1 for clarity. 
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FIG. 11: The best-fit exponent a versus wavevector q for silica of diameter 0.45 /.im measured at 
different s. Points are the fitting results and the lines are guides to the eye. Data are displaced by 
a factor of 1.1 for clarity. 
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FIG. 12: Intensities of 3 mm thickness Gillette shaving foam at sample-to-detector distance s = 103 
mm. The symbols are the data. And the lines are the fittings based on a stretched exponential 
decay I{q) = Cexp(-(^ • qY). 
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FIG. 13: Normalized intermediate scattering functions of 3 mm thickness Gillette shaving foam at 
sample-to-detector distance s = 103 mm and q = 8.18 X 10"^ A. The symbols are the data. And 



the lines are the fittings based on a stretched exponential decay (Eq. 36). 
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FIG. 14: Intensities of Polystyrene 4 /im suspension at different sample-to-detector distances. The 
disks with different colors are the data measured at different sample to detector distance. The solid 
lines are theoretical plots of the transfer function T{q) of the silica suspension for comparison. The 
data are displaced for clearance. 
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FIG. 15: Normalized Autocorrelation Functions of Polystyrene 4 fim suspension measured at dif- 
ferent sample-to-detector distances. 
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FIG. 16: The best-fit decay rate T versus wavevector q for polystyrene of diameter 4 measured 
at different s. Points are tlie fitting results and the lines are guides to the eye. Data are displaced 
for clarity. 
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